$>.4  0  761  54 


7c? 

^  ft-v 


NRL  Memorandum  Report/  4088 


.Electromagnetic  Instabilities  in  a  Focused 

_  Jon  Beam  Propagating  Through  a  ( - 

1 - z-Discharge  Plasma  *  / 

^  - - I 


76)  ;  . ^  . 

~*f  P'  F/>TTINGER’  D'  ^°SHER  ShVKE  AyfcoLDSTEIN  / 


Plasma  Physics  Division 


&L . 

/  Oct 


D  D  C 


NOV  6  1979 


naval  research  laboratory 

Washington,  D,C.  _  , 

Approved  for  public  release;  distribution  unlimited. 


79  U  05 


025 


SECURITY  CLASSIFICATION  of  This  PAGE  (Whan  Data  tnlarad) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  REPORT  NUMBER  2  GOVT  ACCESSION  NO. 

NR  I.  Memorandum  Report  4088 

3  RECIPIENT'S  CATALOG  NUMBER 

i  -  ■  |  ■  ■  - - —  ■ 

4  TITLE  (and  Sublltla) 

KLECTROMAGNFTK  INSTABILITIES  IN  A  FOCUSED 

5  TYPE  OF  REPORT  A  PERIOD  COVERED 

Interium  report  on  a  continuing 
NRL  problem 

ION  BEAM  PROPAGATING  THROUGH  A  z-DISCHARGE 
PLASMA 

6  PERFORMING  ORG  REPORT  NUMBER 

7  author^*; 

P.  F.  Ottinger3,  D.  Mosher,  and  Shyke  A.  Goldstein3 

e  contract  or  grant  numbers; 

9  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Naval  Research  Laboratory  J 

Washington,  D.C.  20375 

10  program  Element  project,  task 

AREA  8  WORK  UNIT  NUMBERS 

FS-77-A-0 1-6021 

NRL  Problem:  67-H02-29B 

11  CONTROLLING  OFFICE  NAME  AND  AODRESS 

12.  REPORT  DATE 

October  15,  1979 

Department  of  Energy 

Washington,  D.C.  20545 

13  NUMBER  OF  PAGES 

25 

A  MONITORING  AGENCY  NAME  A  ADORES Sf//  dlllarant  Irom  Controlling  O/ll ca) 


15.  SECURITY  CLASS,  (ol  thla  raport) 

Unclassified 


15*.  OECL  ASSIFIC  ATION  DOWNGRADING 
SCHEDULE 


17.  DISTRIBUTION  STATEMENT  (ol  tha  abatract  antarad  In  Block  20,  II  dlllarant  Irom  Raport) 


18  SUPPLEMENTARY  NOTES 

This  research  was  supported  by  the  Department  of  Energy  project  number  67-H02-29B. 
a.  Present  address:  JAYCOR,  Alexandria,  Virginia  22304. 


19.  KEY  WORDS  ( Contlnua  on  ravaraa  alda  II  nacaaaary  and  idantlly  by  block  numbar) 

Electromagnetic  Instabilities 
Ion  Beam  Propagation 
Z-Discharge  Plasma  Channel 
Ion  Orbits 


20.  ABSTRACT  (Contlnua  on  ravaraa  alda  II  nacaaasry  and  Idantlly  by  block  numbar) 

^  A  beam-plasma  system  consisting  of  a  focused  light  ion  beam  propagating  through  a 
z-discharge  plasma  is  analyzed  for  electromagnetic  velocity  space  instabilities.  In  particular,  the 
Weibel  (k  ■  D  --0,  It  ■  "*  0)-and  the  Whistler  4L-X  D  "‘*  '0;  *  •  Vj-—  Of  instabilities  are 

studied.  It  is  found  that  unstable  modes  do  not  grow  fast  enough  to  affect  the  beam  propagation 
in  a  pellet-fusion-reactor  chamber.  For  k  •  Vz  >  0,  the  instabilities  convect  with  a  velocity  much 
less  than  the  beam  streaming  velocity^ —  — -. 


DO  ,  r°Z,  1473  EDITION  OF  1  NOV  65  IS  OBSOLETE 

S/N  0  102-0  14-  6601 


1 


SECURITY  CLASSIFICATION  OF  TmiS  PAGE  (Whan  Data  Bntarad) 


CONTENTS 


I.  INTRODUCTION  .  1 

II.  BEAM-PLASMA  EQUILIBRIUM  .  2 

III.  THE  WEIBEL  INSTABILITY  .  4 

A.  Ion  Instability  .  5 

B.  Electron  Instability  .  8 

C.  Summary  for  the  Weibel  Instabilities  .  9 

IV.  THE  WHISTLER  INSTABILITY  .  10 

V.  CONCLUSIONS  .  14 

ACKNOWLEDGMENTS .  16 

APPENDIX  A  .  17 

APPENDIX  B  .  18 

APPENDIX  C  .  20 

REFERENCES  . 22 


ELECTROMAGNETIC  INSTABILITIES  IN  A  FOCUSED 
ION  BEAM  PROPAGATING  THROUGH  A 
Z -DISCHARGE  PLASMA 

I.  INTRODUCTION 

In  an  earlier  paper,1  it  was  shown  that  focused  ion  beams  for  use  in  a  pellet  fusion  device 
can  propagate  axially  down  a  z-discharge  plasma  channel  without  generating  disruptive  micro¬ 
turbulence  due  to  electrostatic  streaming  instabilities.  The  azimuthal  magnetic  field  in  the  z- 
discharge  channel  confines  the  beam  radially  as  it  propagates.  Here  the  analysis  will  be 
extended  to  study  electromagnetic  velocity-space  instabilities.  In  particular,  the  Weibe!  instabil¬ 
ity  (k  B  =  0,  k  ■  V.  =  0)  and  the  Whistler  instability  (k  x  B  =  0,  k  ■  V.  =  0)  are  investi¬ 
gated,  where  k  is  the  wavevector,  B  is  the  azimuthal  magnetic  field  and  V,  =  V.  e.  is  the  axial 
streaming  velocity  of  the  beam. 

The  beam-plasma  system  consists  of  a  focused  ion  beam  (typically  a  5  MeV  proton  beam 
of  50  ns  duration,  0.5  cm  radius,  and  a  current  of  5  x  105/U  propagating  down  the  axis  of  a  z- 
discharge  plasma  channel.2  The  ion  beam  is  focused  at  the  entrance  to  the  plasma  channel  (see 
Figure  1)  with  velocity  components  transverse  to  z  given  by  VJV.  =:  tan  fi  «  I.  A  high 
plasma  density  in  the  channel  (n^  ~  10lscm“3)  insures  good  beam  charge  neutralization.3 
Good  beam  current  neutralization  in  the  interior  of  the  beam  also  occurs,  so  that  the  total  mag¬ 
netic  field  is  comparable  to  that  associated  with  the  preformed  channel  established  before  beam 
injection.  The  beam  current  greatly  exceeds  that  establishing  the  channel  so  the  electron  drift 
velocity  is  approximated  by  Ve  =nhVJnr 
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Hydrodynamic  modeling  of  the  background  plasma4  shows  that  a  uniform  channel  net- 
current  model  is  appropriate  for  the  early  times  associated  with  passage  of  the  beam  front.  This 
is  because  the  low-temperature  channel  is  established  microseconds  before  beam  injection  so 
that  complete  magnetic  diffusion  occurs.  Later  in  the  ion  pulse,  expansion  of  the  beam-heated 
high-temperature  plasma  (T  =  25-50  eV)  reduces  the  magnetic  field  strength  in  the  interior  of 
the  channel.  The  built-up  field  in  the  expanding  cylindrical  shock  wave  is  also  enhanced  by 
significant  current  non-neutralization  in  the  cool  plasma  surrounding  the  beam-heated  channel. 
The  maximum  field  strength  just  outside  the  ion-beam  radius  can  exceed  that  established  by  the 
initial  z-discharge  current  by  a  large  factor.  Thus,  at  late  times  during  beam  passage,  the  mag¬ 
netic  field  distribution  is  closely  approximated  by  a  surface-current  model. 

In  Sec.  II,  equilibrium  models  for  such  a  beam  plasma  system  will  be  described.  In  Secs. 
Ill  and  IV,  the  Wiebel  and  Whistler  instabilities  will  be  investigated  The  conclusions  which 
can  be  drawn  from  this  work  are  summarized  in  Sec.  V. 


II.  BEAM-PLASMA  EQUILIBRIUM 

For  mathematical  convenience,  a  slab  model  will  be  used  for  the  beam-plasma  system. 
This  is  appropriate  for  the  case  at  hand  since  ions  are  injected  into  the  channel  with  small  angu¬ 
lar  momentum  so  that  the  resulting  orbital  motion  occurs  in  a  plane.  At  early  times  in  the 
pulse,  the  net  current  (nearly  equal  to  the  channel  current)  is  uniformly  distributed  across  the 
channel  and  flows  in  the  z  direction.  Thus,  B  -  where 


Bpx/a,  \x\  <  a 


(1) 


Here,  B„  is  the  peak  value  of  the  field  and  "a"  is  the  channel  radius.  If  the  beam  distribution 
function  fh  is  written  as 


n  2  VP 

fh(\  v  v,)  =  —  8(vv)  8(v2  +  v.2 - — 1  +  K) 

it  '  m, 


(2) 
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where  Pt  -  my,  +  eAjc  is  the  axial  canonical  momentum,  A.  is  the  vector  potential  and  nb, 
V.  and  K  are  constants,  then 


nb(x)  - 


i«*.  I*  I  <  rb 

.0.  1*1  >  r„  ’ 


where  rb  -  V ^w0,  K  -  V2  -  Vj,  at2  -  VzwcJa  and  «u(,  -  eBp/m,c.  A  smoothly  falling  density 
profile  may  be  obtained  by  replacing  the  second  delta  function  in  Eq  (2)  by  a  Maxwellian  distri¬ 
bution  function.  The  distribution  function  in  Eq.  (2)  also  states  that  all  beam  ions  cross  the  axis 
at  the  same  angle  and  traverse  the  entire  beam  radius  during  each  betatron  oscillation.  A  more 
complicated  distribution  function  could  be  used  to  model  the  small  spread  in  angles  at  which 
the  ions  cross  the  axis,  however  little  additional  information  is  obtained  for  the  effort.  For 
mathematical  convenience,  the  form  of  the  distribution  function  given  in  Eq  (2)  will  be  used 
here.  It  is  easy  to  show  that  the  fluid  velocity  is  given  by  V  —  Ke.  and  that  in  order  for  the 
beam  to  be  confined  within  the  plasma  channel,  one  must  have  rb  ^  a.  Furthermore,  fb  can 


be  written  in  the  more  convenient  form 


fb(yy,  <£)--—  8(v,,)  8(v^  -  Vl  +  <»2x2),  (4) 

where  vx  =  v0  sin<£,  v,  =  V.  +  v^  cos<l>  and  vj  —  v2  +  (v.  —  V.)2.  Here,  K  is  associated 
with  the  average  streaming  velocity  of  the  beam  ions  and  v0  is  associated  with  the  oscillatory 
betatron  motion  of  the  beam  ions  (v^  <  Ve  -  w0rb).  The  beam-ion  orbit  equations  for  this 
field  geometry  were  solved  in  Ref.  1  and  the  results  are  summarized  in  Appendix  A. 


The  distribution  function  given  in  Eqs.  (2)  or  (4)  provides  an  appropriate  discription  of 
the  ion  beam  at  early  times  in  the  pulse.  At  late  times  Bv  -  0  inside  the  channel  and  the  field 
is  built-up  sharply  at  the  radial  edge  of  the  beam.  Thus  at  late  times  in  the  pulse 

-^-8(v,)  8(vj-  Vl),\x\  <  r„ 
fb(v,v  ve,  </>)  -  q  |x|  >  rb  +  8 
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where  Bv(rb  +  8)  is  sufficiently  large  to  confine  the  beam  and  the  sheath  is  restricted  to  a  thin 
layer  such  that  6  «  rb.  Beam  ions  move  in  straight  line  orbits  inside  the  field-free  channel. 
Within  the  layer  yb  <  |x|  <  yb  +  6,  the  ions  reflect  off  the  magnetic  wall,  reverse  their 
transverse  velocity  and  resume  their  straight  line  trajectories  after  reentering  the  channel.  The 
distribution  function  in  Eq.  (5)  also  results  in  the  uniform  density  profile  of  Eq.  (3)  and  in  a 
fluid  velocity  given  by  V  -  V.e.  inside  the  channel. 

The  background  plasma  provides  complete  charge  and  nearly  complete  current  neutraliza¬ 
tion  of  the  beam.  In  addition,  the  plasma  also  carries  the  z-discharge  current.  The  high  density 
desired  for  good  beam  neutralization  provides  for  a  high  frequency  of  electron-plasma  ion  colli¬ 
sions,  i/,,  shown4  to  be  larger  than  a>ce  inside  the  channel  at  all  times  during  the  pulse.  Thus,  a 
collisional  fluid  model  is  used  for  the  background  plasma  with  the  electrons  drifting  with  velo¬ 
city  \e  =  (nbV.Jnp)  ez. 

III.  THE  WEIBEL  INSTABILITY 

Two  Weibel  instabilities  will  be  investigated,  the  ion  instability  and  the  electron  instabil¬ 
ity,  which  are  respectively  driven  by  the  streaming  of  the  beam  ions  and  the  electron  drift 
motion.  Lee  and  Lampe5  report  for  electron  beams  that  the  Weibel  instability  grows  at  a 
greatly  reduced  rate  when  VJV2  >  where  u>pb  is  the  beam  plasma  frequency  and  is 

the  electron  plasma  frequency.  Molvig6  has  shown  that  is  possible  for  electron-ion  colisions  to 
restore  rapid  growth  of  the  mode.  Although  the  concern  here  is  with  ion  beams  driving  the 
instability,  again  both  beam-thermal  effects  and  collisional  effects  are  important.  Ion  betatron 
motion  will  also  be  important  in  analyzing  the  ion  instability. 
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A.  Ion  Instability 


Consider  first  the  situation  late  in  the  beam  pulse,  with  fb  given  by  Eq.  (5).  For  assumed 
large  rb ,  the  perturbed  distribution  function  with  E(x)  —  E  exp(/Cx), 

e/» 


A- "  ?  J- 


r  0 

_  v  x  B 

»  dT 

E  +  - 

J  _oo 

c 

dv 


exp  i(kvx  -  w )  r. 


(6) 


can  be  integrated  over  t  to  give 

fb\ 


—  ie/m, 

kv. 

r  i  _ i.  r 

/cv, 

1 - 1 

to 

,  «/* 

> 

1 

3 

£,  “T  C.  - 

Oi 

dv, 

*8^7 

(7) 

Here  k  -  kex  and  straight-line  unperturbed  orbits  are  used  since  Bt  -  0  for  |x|  <  rb.  The 
assumption  that  rtb/np  «  1  allows  one  to  write  the  usual  approximate  dispersion  equation,6 
D..  ~  0,  for  the  Weibel  instability  where  |cu|  <  i>,,  and 


D„  ~  0  =  c2k 2  +  y2  +  to2,y/v,  +  /z:. 


(8) 


Here  u  =  iy  for  purely  growing  perturbations,6  v,  —  V,  »  Ve  (electron  drift  motion 

is  ignored  at  present)  and  /..  is  the  beam  contribution  to  D::\ 


L. 


—47 re2  f 

k\ } 

*Ll  4 

a/*' 

m,  J 

> 

1 

V‘  av. 

d\. 


(9) 


Integrating  by  parts  and  then  using  the  calculus  of  residues  to  perform  the  remaining  nontrivial 
integration  results  in 


which  reduces  to 


klV 2 


Ar(Kg  +  2F.) 

(1  +  y2/k2V2)'12 


kv}/  y0 

(1  +  y2/k2V%)112 


7k7V2Vp 

y 


(10) 


2  F  F,2y 

— -  +  — — 
**  kVl 

,  y  «  kV p 

k2V2 

1  , 

, 

y  »  kV0 

y 

(11) 
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for  V,  >  Vp  and  the  limits  shown.  Solving  Eq.  (8)  for  the  growth  rate,  y,  one  obtains 


<k2y?*  0  <k<k, 

_  I  K- 1  |  2w*tv*vl\ 

I  -}v,  I  (  -S*i  )'*'<*<*•• 

2<opbc(2VJV^lHk0-  k),  k  =  k0 


kp  -  («>,  y?/2il2ajpVp);  k0  -  (b>pb/c)  (2  VJ  Vp) 1 


The  peak  growth  rate  is  given  by 


(  2  K2 


I  ^ 


at  k  —  k„ 


For  V",  ~  10  3  and  V  K  >  10  \  y‘p  -  106  sec  1  at  late  times  in  the  beam  pulse  (note  that 
yp  is  actually  overestimated  here  since  kprb  «  1).  Thus,  no  significant  growth  can  occur  since 
yPTb  >  0.05;  here  Th  is  the  beam  pulse  length. 


At  early  times  in  the  beam  pulse,  v,  is  larger  due  to  low  channel  temperatures  so  that  one 
might  expect  the  growth  rate  to  be  larger.  It  will  be  found,  however,  that  by  including  the 
betatron  motion  of  beam  ions  in  the  analysis,  the  perturbation  is  stabilized.  In  this  case,  fb  is 
given  by  Eq.  (4),  the  ion  orbits  are  found  in  Appendix  A,  and  /M  is  given  by 


-2e  dfb 


\e  oJh  r" 

7  ^7  J_„  dT  Vx'Ex(x')  +  (v:  -  K)  E.(x') 


iV:vx'  dE; 

a*  dx' 


exp  (/tor) 


(14) 


It  is  now  assumed  that  E,(x)  =  E.  cos  kx  where  k  is  restricted  to  a  discrete  set  of  values 
by  boundary  conditions  at  rh  with  krh  >  1.  This  choice  of  £.(*)  is  reasonble  since  np  and  nh 
are  both  uniform  and  nh  «  np.  Thus,  the  mode  is  expected  to  closely  resemble  the  back¬ 
ground  plasma  eigenmode  which  is  sinusoidal.  Here  again,  the  background  plasma  is  treated  as 
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a  collisional  fluid  since  ve  ><uiv.  An  approximate  algebraic  dispersion  equation  is  obtained  by 
taking  a  weighted  spatial  average  of  DnE.(x)  where  D::  is  a  diffe'ential  operator  in  x  defined  by 

D::  *  -c2-^-j-  +  y2  +  ~~  -  f  vjb\d\ 

dxz  v,  J 

and  ft, |  is  given  in  Eq.  (14)  (£,  =  0).  The  dispersion  equation  is  then 


0  =  D„  =  — — -  f  dx  D.,E.(x)  cos  kx 
rhE: 


(15) 


Using  the  orbit  equations  for  x'  and  v’  found  in  Appendix  A  and  the  Bessel  function  iden¬ 


tities 


exp(±  iz  sin  0)  -  J0(z )  +  2  Y.J 2,(z)  cos  219 


i- 1 


and 


±  2/  ^  y2/+I (2 )  sin  (21  +  1)  9, 

1- 0 


exp(/z  cos  9 )  =»  J0(z)  +  2  £  f  J,(z)  cos  19, 


the  t  integration  in  Eq.  (14)  can  be  performed.  Thus,  Eq.  (15)  takes  the  form 


(16) 


(17) 


0  —  c2k 2  +  y2  +  — — -  +  /„ 

v, 

where  in  terms  of  aqU.  \x,  v.)  and  a2(x,  \x,  v.) 


(18) 


7  Cr«  ,  ,  r  ,  1 

-  -  I  dx  cos  kx  I  d  v  v. - - 

nb  J-x  J  8v  j 


kV. 


a  1  +  a  • 


(19) 


The  expressions  for  a,  and  a2  are  given  in  Appendix  B.  Terms  involving  resonances  at  higher 
harmonics  of  a>0  [  i.e.  (<o  -  mw„)_l  for  m  ^  0]  are  not  considered.  The  term  involving  a,  in 
Eq.  (19)  is  the  source  of  instability.  Here,  however,  because  of  the  betatron  motion,  this  term 
vanishes  and  the  mode  is  stable  at  early  times  in  the  pulse  (see  Appendix  B).  This  same  result 
is  obtained  for  odd  eigenfunctions,  E.ix)  =  E,  sin  fcx.  Mathematically  the  term  vanishes 
because  the  integrand  involving  a!  is  an  odd  function  of  x  (see  Appendix  B).  Physically  no 
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radial  current  bunching  can  occur  since  each  ion  travels  radially  back  and  forth  across  the  beam 
profile  as  it  follows  its  betatron  orbit.  However,  azimuthal  current  bunching  can  occur  as  will 
be  discussed  in  Sec.  IV. 


B.  Electron  Instability 


The  only  other  source  of  energy  available  to  drive  the  Weibel  instability  at  early  times  in 
the  pulse  is  the  drifting  electron  background.  The  drift  velocity,  however,  is  actually  subther- 
mal  since  even  before  beam  heating  occurs  Ve  —  VJ np  <  ue  —  (7’/m(.)l/2.  Using  a  warm 
collisional  model  (v,,  >  wie  at  all  times),  the  dispersion  equation  for  the  electron-Weibel  insta¬ 
bility  is  derived  in  Appendix  C.  Setting  y  =  0  in  Eq.  (C7),  it  is  found  that  y  >  0  for 
0  <  k  <  k„-  ojpiVe/y/2u:c ,  where  u,  «  ( T/m,)u 2.  Thus  for  k  <  k„  and  w2p,  >  yv„  Eq.  (C8) 
reduces  to 


0 


c2k2  + 


<^n,y 

V, 


wplk2y2(y2  +  v,y  +  k2u,2)  _ 

(y4  +  v,y 3  +  3k2«,2y1  +  2 v,k2u}y  +  2 k4u,4)  ' 


where  the  ion  beam  contribution  to  Eq.  (20)  is  ignorable. 


(20) 


At  early  times  in  the  pulse,  wp,ujcv,  <  1  and 


y 


(k2V2vy\  0  <  k  <  k\ 

lay  VJe,  *i  <  k  <  k0 


J2  v,cVe 

a>,,u, 


(k0  —  k),  k  ~  k0 


(21) 


where  k,  -  (a> Jc)  (a.„V'(,/ci/,),/2.  Here  T  =  4  eV,  n„  =  2  x  1018  cm  3  and  nb/np  =  SxlO  4, 
so  that  the  peak  growth  rate  is  on  the  order  of  yr  -  a >PY Jc  -  4.2  rb\ 


At  later  times  in  the  pulse,  the  beam  has  heated  the  plasma;  then  wpp,/cv,  >  1  and 


y 


(k2V2v)  l/3, 

V2vJ2u2 


kVe  lk2c2u2 

2  1  _  a)  j,  V 2 


0  <  k  <  k, 
k|  <  k  <  k2 

1/2 

,  k2  <  k  ^  k„ 


(22) 
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where  k,  -  v,  Vt2/u,3,  k2-v,Ve/u,2  and  Ve  <  u,.  Now  y/  -  w„F,2/(23/3cu,)  =  I.Ot*"*  for 
T  —  38  eF,  np  —  l  x  1018  cm  3  and  nb/np  —  1  x  10~3.  Thus,  yeptb  decreases  by  a  factor  of 
four  as  the  beam  passes  through  and  heats  the  plasma.  In  fact,  <op,u,/ci>,  -  1  at  r  ^  10.0  ns 
into  the  pulse  so  that  a  total  of  about  1.6  e-folds  occur  during  the  transit  of  the  beam. 


Note  that  t„  is  the  appropriate  time  scale  for  instability  growth  since  Ve  -  0  except  dur¬ 
ing  the  passage  of  the  beam.  Furthermore,  for  small  k,  (i.e.  k?  <  2wjJ3c J),  it  can  be  shown 
from  Eq.  (C9)  that 


dur 
~dk ~ 


3k?  V? 


2y„2 


-  <  k  «  y. 


and 


:  yn 


k?  v} 

2y0 


(23) 


(24) 


where  y„  -=  tup,VeU..  Thus,  the  unstable  mode  convects  axially,  but  with  a  group  velocity,  vs, 
which  is  much  slower  than  V.  or  Ve.  At  any  given  point  in  the  plasma  the  mode  grows  only  for 
a  time  of  order  rb.  For  larger  k.  and  fixed  kx  the  mode  transforms  into  the  electrostatic 
streaming  instability  which  was  found  to  be  stable  in  Ref.  1. 


C.  Summary  for  the  Weibel  Instabilities 


In  summary,  it  is  found  that  the  betatron  motion  of  the  beam  ions  stabilizes  the  ion- 
Weibel  instability  at  the  beam  front  while  growth  is  too  slow  at  the  tail  of  the  beam  to  allow  for 
even  one  e-fold  (yprb  >  0.05).  The  electron-Weiber  instability,  on  the  otherhand,  grows 
fastest  (ypTb  =  4.2)  at  the  front  of  the  beam  where  the  plasma  is  relatively  cold.  At  the  tail  of 
the  beam  y/r*  =  1.0.  Although  the  eiectron-Weibel  instability  grows  faster  than  the  ion- 
Weibel  instability,  it  also  is  not  expected  to  grow  to  a  level  which  could  seriously  affect  beam 


propagation. 
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IV.  THE  WHISTLER  INSTABILITY 


The  Whistler  mode  (k  x  B  =s  0,  k  ■  V,  —  0)  like  the  Weibel  instability  can  be  driven 
unstable  by  particle  streaming.  However,  the  wavevector  k  —  kyey  (i.e.  e0  in  cylindrical 
geometry)  is  perpendicular  to  the  direction  of  the  betatron  motion  of  the  beam  ions.  Hence, 
unlike  the  Weibel  instability,  the  Whistler  instability  cannot  be  stabilized  by  the  beam  ion  beta¬ 
tron  motion.  However,  a  small  spread  in  v„  (angular  momemtum)  can  reduce  the  growth  rate 
significantly.  The  electron  drift  velocity  is  ignored  when  considering  the  ion-Whistler  instability 
since  V,  »  Ve. 


Consider  first  the  situation  late  in  the  beam  pulse  with  fh  given  by  Eq.  (5).  If  it  is 
assumed  for  the  moment  that  rh  is  very  large,  then  the  perturbed  distribution  function  is  given 
by 


^  -te  „  ,  v  x  B  d/h  „  ,  , 

Jb\  ~  -  E  +  -  '  —  /(«a  ~  k„vv).  (25) 

m.  c 

Using  this  expression  for  fhh  the  perturbed  current  is  easily  calculated  and  used  to  derive  the 


dispersion  equation. 


„  ,  ,  ,  oij/y  ,(  k2(V?  +  Vlf 2) 

D::  «  c2k2y  +  y2  + - +  (ojjl - - — — -  =0.  (26) 

v,  'I  y2 

Again  nh/nn  «  1  was  used  in  deriving  Eq.  (26)  where  a>  =  iy  for  purely  growing  perturba¬ 
tions  and  Ve  «  V.  (electron  drift  motion  is  ignored  at  present). 


Solving  Eq.  (26)  yields 

*y V2.  0  <  ky<  k{ 

y  =  (k2V2u2hv,/oi2)'n,  k\  <  ky  <  k2  ,  (27) 

‘ OpbV-Jc ,  ky  >  k  2 

where 

K  >  Vp  k\  =  o)2bv,lul>2Vz\  k2  =  (u>Jc)  (ionhV,/v,c)U2. 
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The  peak  growth  rate  (for  kt  >  k2  =  25.0  cm1)  is  given  by  y'  -  <ophVz/c  =  2  x  102  for 
a  beam  with  no  spread  in  v,. 


This  result  is  only  slightly  modified  when  finite  geometry  effects  and  the  betatron  motion 
of  the  beam  ions  and  included.  Proceeding  as  in  Appendix  B  with  k  -  kxix  +  kyit  and  fb 
given  in  Eq.  (4),  one  obtains 


D::  -  c2(k2  +  fc2)  +y2+^+  <oi 


pb 


HKO  -  Mat 


(28) 


where  F(kxr„)  is  defined  in  Eq.  (B6).  The  quantity  G(kxrb)  is  defined  by 


c s  cos  x\j„(X)  J0HZ) 


+  2  Z  h n(X)  JflZ) 


,  dX 


(29) 


JZ-2m 

with  X  “  kxx,  Z(X)  =  (kxrb/2)  (1  —  X2/kxrb)]  '  and  N  =  rb  +  (sin  2kxrb)/2kx.  Note  that 
G(kx  =  0)  =  1.  fand  C  are  only  geometrical  factors  and  do  not  modify  the  structure  of  the 
dispersion  equation.  For  large  ky  the  peak  growth  rate  now  becomes 
y'n  =  (a>PbVzlc)  I G(kxrb))xn,  which,  aside  from  the  geometrical  correction,  is  identical  to  y 
found  in  Eq.  (27). 


If  the  beam  has  a  small  spread  in  v^,  the  beam  can  be  modeled  by  the  distribution  func¬ 
tion 

fb  “  ~t”V  S(v»  “  exp(-v2/K2),  (30) 

where  Vy  is  the  thermal  velocity  of  the  beam  ions.  Since  the  y  motion  is  unaffected  by  Bxx  the 
distribution  in  \y  will  be  the  same  at  all  points  inside  the  channel.  For  convenience  finite 
geometry  effects  and  the  effects  of  the  betatron  motion  of  the  beam  ions  are  ignored  here  since 
they  have  little  affect  on  the  Whistler  mode.  Substituting  the  distribution  function  of  Eq.  (30) 
into  Eq.  (25)  yields  an  expression  for  the  perturbed  distribution  function,  fb].  The  perturbed 
current  jb ,  can  then  be  calculated  and  the  dispersion  equation  derived. 


II 


In  Eq.  (31),  oj  =  iy  and 
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D..  —  c2k2  +  y2  +  oj2y/v,  +  /..  —  0, 


io:b  r°°  |  2Ar,vv(  V?  +  Vn/2)  I  ,  . 

4  *  -zr,  /-  *>1' +  -*;»,)  |  “p<_v-/  ^  ’  021 

Here,  the  \e  and  d>  integrations  are  trivial  and  have  already  been  carried  out.  The  quantity  /„ 
is  easily  written  in  terms  of  the  usual  plasma  dispersion  function,  Z(£) 8.  However,  here  £  is 
pure  imaginary  so  that 


,  ,  .  2K’l.  nS» 

'=  -  i  -  -pr|'  - 1±  «»  * 


where  V.  >  V B  and  erf(x)  is  the  usual  error  function.9  For  y  >  kV^  Eq.  (31)  reduces  to  Eq. 
(26),  but  for  y  <  AcFv  the  dispersion  equation  becomes 


0  =  c2k2  +  yJ+^  +  uil- 


reduced  if  V2  >  2v ,(orbcV:/a>pi2.  Since  this  is  the  case  here 


^ Y  ,  2y2 

/* \A\ 

kVy  k2v2  ] 

<inb  Vjc,  for  the  case  with 

Vv  —  0  is 

kV:, 

y~  (k2V:2p,(o2b/to2)u\ 

2u?bp,  V?  _  c2k2  V2 
<o2p,V2  2<tipt,  V2 


0  <  Ac  <  Ac, 

Ac ,  <  k  <  k2 

k-,  <  k  s?  k„ 


where  Ac,  =  w2bvJa>2Vz,  k2  =  Ac,  V?/V 2  and  k„  «  V2  o>ph  VJcVy. 

The  peak  growth  rate  is  now 

Yn  “  (wnh  V-Jc)  (2v,a>pbcV2l  V2w2,)  <  oj pbVjc, 

where  y 'p  varies  in  time  as  T~i/2(t)  through  v,.  The  number  of  e-folds,  8  during  the  beam 


transit  at  any  point  in  the  plasma  is  given  by 

8  -  f’*  y-O)  d,  -  y;(0) 


dt 

(1  +  bTt/rb)il2' 
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where  it  is  demonstrated  later  than  no  significant  wave  convection  occurs.  Here  A  7" 
-  [ 7'(t<>)  -  HO)]/ HO)  for  T in  eFand  >;( 0)  is  the  peak  growth  rate  at  the  front  of  the  beam. 
A  linear  rise  in  temperature,  T(t)/T(0)  —  1  +  A  Tt/rb,  agrees  well  with  results  of  previous 
work.3  Thus 


2y»(0)  Tb 


(1  +  AD1 


For  T(rb)  40  eV,  7X0)  -  4  eV,  rb  -  50  ns,  VJV.  >  0.08  and  nb/nr  -  10"3,  8  <  1.0  e-folds. 
A  spread  in  Vv  then  reduces  the  growth  of  the  ion-Whistler  instability  to  a  tolerable  level. 


Finally  it  can  also  be  shown  that  the  ion-Whistler  instability  does  not  convect  with  the 
beam  when  k,  >  0.  This  needs  to  be  verified  in  order  to  justify  using  rb  as  the  appropriate 
time  scale.  Taking  Eq.  (5)  for  fb  and  setting  k  -  *,e,  +  k.e.  this  dispersion  equation  becomes 


0  —  c2k2  —  (to,  +  iy )•’  f  — —  (to,  +  ly)  + 

v , 

where  to  =  to,  +  ry  and 


(38) 


to3,  to2  +  k2(y?  +  2  V.  Fflcostt  +  Fj  cos2<<>) 

/,.  —  — —  I  d<b - — — - — — - — : -  (39) 

27 r  (to  -  k.V.  -  k:V0cos<t>)2 

Here  the  spread  in  v,  is  neglected  since  it  will  have  iiti'e  effect  on  the  axial  group  velocity  of 
the  perturbation.  For  mathematical  convenience  finite  geometry  effects  and  the  betatron 
motion  of  the  beam  ions  are  also  ignored.  The  results  of  the  present  calculation  is  actually  an 
upper  limit  on  the  group  velocity  since  the  betatron  motion  of  the  beam  ions  will  tend  to  wash 
out  any  disturbance  moving  axially  on  the  beam.  For  |ai,  —  k.V,\  >  k.V B  the  denominator  of 
the  integrand  in  Eq.  (39)  can  be  expanded  and  the  integration  is  trivial.  However,  if 
|eu,  -  /c 2 F, |  <  k.VB  the  integration  is  most  easily  done  using  the  calculus  of  residues.  For 


small  k.  and  V0  <  V.,  |«a,  -  k.V. \  >  k.V0  is  appropriate  for  the  ion-Whistler  instability.  Then 
integrating  Eq.  (39),  Eq.  (38)  can  be  solved  for  wr  and  y  yielding 
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y  -  (  v}  +  V\l 2) 1/2  -  A:,  V,.  (40) 

-  y  k?V2dtopb.  (41) 

Clearly,  \K  —  do>r/dk.  <  V.  for  /c,  <  upb/3c.  Thus  the  ion-Whistler  instability  does  convect  in 
the  axial  direction  for  k.  <  wpb/3c ,  but  with  a  group  velocity  slower  than  the  beam  streaming 
velocity.  For  k.  >  upb/3c  the  calculus  of  residues  can  be  used  to  evaluate  Eq.  (39),  however, 
the  mode  is  then  basically  an  electrostatic  two  stream  mode.  This  mode  has  already  been 
shown  to  be  stable  in  Ref.  1. 

In  summary  it  is  found  that,  as  expected,  the  betatron  motion  of  the  beam  ions  does  not 
affect  the  ion-Whistler  instability.  The  peak  growth  rate,  however,  can  be  reduced  to  a  toler¬ 
able  level  by  the  presence  of  a  spread  in  v^  (angular  momentum  in  cylindrical  geometry).  The 
spread  in  v,,  known  to  be  present  at  injection  in  typical  experiments2  is  sufficient  to  reduce  the 
number  of  e-folds  to  less  than  1.0.  Furthermore,  the  mode  convects  axially  at  a  group  velocity 
less  than  V..  The  electron-Whistler  instability  was  not  considered,  since  it  will  have  properties 
similar  to  the  properties  of  the  electron-Weibel  instability  (discussed  in  Section  II. B)  for  such  a 
highly  collisional  plasma  ( ve  >  wce). 

V.  CONCLUSIONS 

The  purpose  of  this  paper  was  to  study  electromagnetic  velocity-space  instabilities  gen¬ 
erated  by  a  focused  ion  beam  propagating  through  a  z-discharge  plasma.  In  particular,  the 
Weibel  instability  (k  •  B  =  0,  k  •  V2  =e  0)  and  the  Whistler  instability  (k  x  B  =  0,  k  •  V;  ~ 
0)  were  investigated.  This  work  is  an  extension  of  the  work  in  Ref.  1 ,  where  electrostatic  insta¬ 
bilities  were  investigated. 

The  ion-Weibel  instability  (driven  by  the  streaming  energy  of  the  beam  ions)  is  found  to 
be  stabilized  by  the  betatron  motion  of  the  beam  ions  at  the  front  of  the  beam.  At  the  tail  of 
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the  beam,  where  beam  ions  follow  straight  line  orbits,  the  growth  of  the  ion-Weibel  instability 
is  two  slow  to  allow  for  even  one  e-fold  (y* rb  =  0.05)  during  the  transit  of  the  beam.  The 
electron-Weibel  instability  (driven  by  the  drifting  plasma  electrons),  on  the  other  hand,  grows 
fastest  (y/rA  =  4.2)  at  the  front  of  the  beam  where  the  plasma  is  relatively  cold  (T=  4  eV). 
At  the  tail  of  the  beam,  where  T  rises  to  about  25-50  eV,  y'Tb  =  1.0.  Although  the  electron- 
Weibel  instability  grows  faster  than  the  ion-Weibel  instability,  it  also  is  not  expected  to  grow  to 
a  level  which  could  drastically  affect  beam  propagation.  Only  1.6  e-folds  will  occur  during  the 
transit  of  the  beam.  It  has  also  been  shown  for  0  <  k?  <  2o»(),2/3c2,  where  k  -  kxix  +  k2e„ 
that  the  electron-Weibel  instability  does  convect  axially  but  at  a  group  velocity  much  less  than 
the  beam  velocity  (v„  —  3fc/c2  F,/2<»>2,  <  V,  «  V2).  Thus  the  appropriate  growth  period  for 
calculating  the  number  of  e-folds  is  just  the  beam  transit  time,  rb ,  and  growth  occurs  mostly  at 
the  tail  of  the  beam. 

Because  the  plasma  is  highly  collisional  (vc>  wa,  at  all  times),  the  electron-Whistler  ins¬ 
tability  will  have  properties  similar  to  those  of  the  electron-Weibel  instability.  Thus  it  is  also 
not  expected  to  grow  to  a  level  which  could  drastically  affect  beam  propagation. 

The  ion-Whistler  instability,  as  expected,  is  not  stabilized  by  the  betatron  motion  of  the 
beam  ions.  The  peak  growth  rate,  however,  can  be  reduced  to  a  tolerable  level  by  the  presence 
of  a  spread  in  v,  (angular  momentum  in  cylindrical  geometry).  Furthermore,  for  k.  <  gj^/3 c 
the  mode  also  convects  axially  at  a  group  velocity  less  than  K(v„  -  3k.cV./<n nh) .  Thus  rh  is 
again  the  correct  time  scale  and  only  1.0  e-fold  are  expected  to  occur  during  the  transit  of  the 
beam  in  typical  experiments. 

From  these  results  and  from  the  results  for  electrostatic  instabilities  in  Ref.  1,  it  can  be 
concluded  that  it  is  possible  to  propagate  a  focused  ion  beam,  appropriate  for  a  pellet  fusion 
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device,  through  a  z-discharge  plasma  channel  without  generating  significant  growth  of  microin¬ 
stabilities. 
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APPENDIX  A 


For  the  magnetic  field  configuration  given  in  Eq.  (1)  for  \x\  <  a  the  beam  ion  orbit  equa¬ 
tions  were  solved  in  paper  I.  The  resulting  ion  orbits  are  given  below: 


and 


x'  —  x  cos  u>„t  +  —  sin  u)0t, 


y  -  y  +  v„r 


8F, 


z'  -  z  + 
2 


«0 

v2 

4  K 

—  x‘ 


r2- 


(On 


XV, 


sin  2<d„t  —  — —  (cos  2o>0t  —  1), 

r  , 


(Al) 

(A2) 

(A3) 


vx'  -  vx  cos  o)„  t  —  xw0  sin  to0r, 

v>' 
v; 


V,  -  V,  - 


4K 


'y  -  V 

,2 


.2 


(COS  2(0 0T  —  1) 


(O0XVx 

+  — tt: —  sin  2<o0t, 


2  V, 

where  <o0  and  Vz  are  the  same  as  defined  in  Section  II. 


(A4) 

(A5) 


(A6) 


17 


APPENDIX  B 


The  expression  for  ct)  and  a2  found  in  Eq.  (19)  are 


«i 


1 


v,W) 


N 

-xa>J{(Y) 


j  j„(y)  +  y2(n 


j  J„(X)  +  J2(X) 


“2“  ^ 


■f  •/„(*)  y0(T)  +  2  i  (-d* /,„(*)  y2„(y) 


n  — 1 


8K 


~ t  -  *2  km  y„(  y)  -  VAX)  j j(  y> 


+  270(K)  y2(AO  +  4  £  (-1)" 

/i  — I 


(Bl) 


Ju(x)  i2,+,(y)  +  -  y2„(y)  r2„+lm 


i!^r1  y,(x)  /,(n  +  i  (-n|/2„+1u>  jwy)  +  j2„_,u)  y2„+1(y>]JJ.  (B2) 


2K 

where  X  -  kx.  Y  -  <rv,/o>„  and  )V  -  r*  +  sin  Ikrjlk. 


In  order  to  show  that 


/ 1  «■ - —  (  *  dx  cos  kx  f  d\ 

nh  J~rb  * J 


kV:y2  d/t 


y  dvj 


“i  “  °- 


(B3) 


first  write  \x  —  vg  sin  <t  and  v.  —  V:  +  vg  cos  d>  and  then  perform  the  v>  and  <b  integration 
yielding 


/,  - 


-»j>kV,2 


yN 

x  Ly,m 


f_*  dx  cos  kx  J  dvg 


88(v  j  -  y|  4-  ft>2x2) 


T  y-i/2(z>  +  Jm(Z) 


6vj 


Jm(Z) 


—  XtO0J  2/2  (Z) 


jjjx)  +  jax) 


(B4) 
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where  Z  —  k\J2a>0.  Integrating  by  parts  in  to  remove  the  derivative  of  the  delta  function 
and  then  performing  the  v@  integration  results  in 


r*'»  J  •/>(*) 

-  — —  I  dX  cos  A1 — - — 
AyN a)0  J~krb  I  Z 


f  y-t/ifz)  +  jm(z) 


JmlZ) 


+  Jx{X) 


y  y'-1/2(z)  y,/,(z)  +  y  y.l/2(z)  y'1/2<z> 


+  J'yz^Z)  J i/2(Z)  +  J'm(Z) 


-  4  ywz) 


1/2* 


y  y#(z)  +  y2(*) 


Z-2(x) 


where  2(x )  -  (krb/2)  (1  -  x2/k2rb)in.  Since  the  integrand  is  an  odd  function  of  x,  /, 


(BS) 

0. 


The  integral 

/2  - - —  J*  ”  dx  cos  fcx  f  </3v 

is  more  complicated  and  the  associated  term  in  the  dispersion  relation  does  not  contribute  to 
instability;  hence  it  will  suffice  to  write  /2  in  terms  of  the  function  F(krb)  as  defined  in  (B6). 


3A|  = 

!  2vfl 


v,a2-rfr  =  up»  F(krJ 


(B6) 


APPENDIX  C 


For  a  warm  collisional  fluid,  the  continuity  equation  and  momentum  transport  equations 
(Te  =  T,  and  v,  —  mtvtlm,)  are 


Hr -  v  •(«-’.)-  o. 


dV“  .  Qa 
-dT+v“Vy“-^: 


E  + 


»„  x  B 


(Cl) 


- —  Vna  -  va(\a  -  \p),  (C2) 

^  a 

for  a  -  i,e.  Linearizing  Eqs.  (Cl)  and  (C2)  for  perturbations  with  exp  i(kx  -  wt)  dependence 
and  solving  for  the  perturbed  current  results  in 


where 


\e+l 


I  •  E, 

47 r 


(C3) 


1  + 


ft,2 


/ft] 


-  4c  -  -/fcK,n,/ yftj2. 

/  _L  _  *2^2»3 

y2ft|  ’ 

where  to  —  /'y,  fl|  —  y  +  k2u2/y ,  ft 2  “  y*',  +  fc2u,2  +  li>,k2u2/y  +  k4u*/y2, 
n3  —  y  +  i>,  +  k2u2/y.  In  deriving  expressions  (C4)-(C6)  it  was  assumed  that  y  <  ,i/r. 


(C4) 

(C5) 

(C6) 

and 


Here  the  electron  streaming  is  driving  the  instability  and  nf  =  n,  (unlike  njnc  «  1), 
thus  the  complete  dispersion  equation, 

DxxD!:  -  Dx:D:x  -  0,  (C7) 

must  be  used.  This  results  in  the  following  dispersion  equation 
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0  - 


,1  +  =A,  +  £l' 
n3  n,2 


>2  +  c2k:+  y^k  _ 


yil: 


n,< 


(C8) 


In  order  to  study  the  axial  convection  of  the  electron  Weibel  instability,  the  dispersion 
equation  must  be  rederived  with  >  0.  For  mathematical  convenience  thermal  effects  are 
ignored  (7  —  0  in  Eq.  (C2)).  In  this  case,  for  large  kx ,  the  dispersion  equation  is 


0  —  Ic2*,2  —  <u2  — 


'“jV 


c2k  2  _  ik*  +  iv) 


(DOJ  V , 


,  ,  ,  , 

k:kxc  + - 


V , 


kzkxc  + 


•ufrxK 

V , 


1  - 


ik2  V ev , 


(C9) 


where  at1  —  at  —  kz  Ve  and  now  w  —  wr  +  iy. 
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Figure  1  -  Typical  ion  trajectory  in  the  confining  azimuthal  magnetic  field  of  the  z-discharge  channel  after  entering  the 
channel  from  the  focussing  region  on  the  left.  Here  Fq  is  the  speed  of  the  beam  ion.  a  in  the  injection  angle,  r,  is  the 
spot  size  at  injection  and  a  is  the  channel  radius 


